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A variant on the kappa opioid receptor gene (OPRK1) is 
associated with stress response and related drug 
craving, limbic brain activation and cocaine relapse 
risk 

K Xu\ D Seo\ C Hodgkinson 2 , Y Hu 3 , D Goldman 2 and R Sinha 1 ' 4 ' 5 

Stress increases drug craving and relapse risk. The kappa opioid receptor gene (OPRK1) mediates stress responses. Here, we 
examined whether the OPRK1 rs6989250 C>G affects stress-induced cocaine craving and Cortisol responses, subsequent 
cocaine relapse risk and the neural response to stress using functional magnetic resonance imaging (fMRI) in cocaine 
dependence. Sixty-seven treatment-engaged, abstinent cocaine-dependent African-Americans were genotyped (CG: N= 10; CC: 
N= 57) and participated in a 3-day experiment in which they were exposed to personalized script-driven imagery of stress, drug 
cues and neutral scenarios, one condition per day, randomly assigned and counterbalanced across subjects. Repeated 
measures of craving and Cortisol were obtained. The subjects were followed prospectively for 90 days to assess relapse risk. A 
follow-up preliminary fMRI experiment assessed neural responses to stress, drug cue and neutral conditions in matched CG 
(A/=5) and CC (A/=8) subgroups. We found greater stress-induced craving (P= 0.019), higher Cortisol during stress and cue 
relative to the neutral condition (P's< 0.003), and increased cocaine relapse risk (P= 0.0075) in the CG compared with the CC 
group. The CG relative to the CC group also showed greater activation of limbic and midbrain regions during stress and cues 
relative to the neutral condition with additional stress-induced activation in the right amygdala/hippocampus (P<0.05, whole- 
brain corrected). These results suggest that OPRK1 is associated with stress-induced craving and Cortisol, hyperactive 
hypothalamus/thalamus-midbrain-cerebellum responses, and also associated with greater subsequent cocaine relapse risk. 
Future studies to replicate these findings in a larger sample size are warranted. 
Translational Psychiatry (2013) 3, e292; doi:1 0.1 038/tp.201 3.62; published online 20 August 2013 



Introduction 

Relapse to drug use is the central problem in the treatment of 
addiction. 1 It is well known that acute and chronic stress 
increases drug craving and also contributes to high rates of 
addiction relapse. 2,3 However, the remarkable individual 
variation observed in stress and drug craving suggests that 
genetic variants may contribute to the differences in response 
to stress, craving and relapse risk. 

Converging evidence shows that the kappa opioid receptor 
(KOR) has an important role in mediating the stress response 
and in stress-related addictive processes. 4-8 Activation of 
the KOR through the administration of dynorphin, a KOR 
endogenous neuropeptide, produces stress aversive symp- 
toms in rats. 6 ' 7 ' 9 ' 10 Physical stress and social defeat elevate 
dynorphin levels, 11-14 suggesting that activation of the KOR 
system increases both the physical and psychological stress 
responses. KOR activation also provokes stress-induced drug 
reward. 15 The administration of a KOR agonist mimics the 
forced swim stress-induced potentiation of cocaine place 
preference. 16 In healthy human subjects, the administration of 



a KOR agonist produces dysphoria and anhedonia, 5 ' 6 which 
resemble withdrawal 'crash' symptoms in cocaine-dependent 
patients. Clinically, stress and dysphoria are closely asso- 
ciated with craving and relapse risk. 2 ' 17-20 Genetic variants 
of the prodynorphin gene are associated with negative 
emotions and alcohol dependence 21 Thus, genetic variants 
of OPRK1 may modulate individual stress responses, stress- 
induced craving and the relapse risk associated with cocaine 
dependence. 

KOR has significant effects on the brain regions involved in 
addiction. 15,22 Repeated administration of a KOR agonist can 
increase dopamine release and cocaine-induced locomotor 
activity, akin to cross-sensitization of stress and drug 
reward 23 Dynorphin mRNA expression and KOR binding 
are elevated in dopamine-rich regions in the post-mortem 
brain tissue of cocaine-dependent patients 24 These findings 
suggest that the activation of KOR produces drug reward 
effects by increasing dopamine release. Thus, it is possible 
that KOR may affect stress-related addictive phenotypes, 
such as craving and relapse. 
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Genetic variants of OPRK1 have previously been linked to 
substance dependence with inconsistent results. 8 ' 21 ' 25-28 
One single-nucleotide polymorphism (SNP) in intron 2 of 
OPRK1 has been associated with alcohol dependence, 25 
whereas three SNPs in the coding region show no association 
with this condition. 29 The frequency of OPRK1 haplotypes 
with intronic SNPs was found to be higher in alcohol- 
dependent subjects compared with controls. 26 A recent study 
found that the non-coding SNP rs99791 7 was associated with 
alcohol-induced sedation. 30 Furthermore, an SNP in a non- 
coding region was associated with opioid dependence in a 
Caucasian population. 27 These data suggest that the reg- 
ulatory region of OPRK1 may be important for substance 
dependence. In support of this notion, a promoter SNP was 
identified in alcohol-avoiding rats, 31 and a promoter 800-bp 
indel variant of OPRK1 was found to affect gene expression 
in vitro. 32 The longer allele showed 50% lower expression 
than the shorter allele. Interestingly, this functional promoter 
variant has been associated with alcohol dependence in 
families. 32 However, no previous research has assessed the 
effects of OPRK1 variants and stress, drug craving and 
relapse risk in a cocaine-dependent population. 

In this study, we examined whether the SNP rs6989250 
C>G on the 5'-regulatory region of OPRK1 affected stress- 
induced cocaine craving and Cortisol in the laboratory and 
subsequent cocaine relapse risk. We focused on this 
particular SNP because: (1) it is located in the 5'-regulatory 
region of OPRK1, and is close to the functional indel variant 
described above and (2) the rare G allele creates an NF-kB 
transcription binding site, suggesting a potential alteration of 
OPRK1 function. Because G allele at rs6989250 is rare 
in Caucasians but abundant in Africans (http://hapmap. 
ncbi.nlm.nih.gov/cgi-perl/gbrowse/hapmap28_B36/#search), 
we tested the association of rs6989250 with stress, drug 
craving and relapse risk in an African-American sample of 
cocaine-dependent patients. As greater stress reactivity 
would also suggest higher neural activity in stress-sensitive 
limbic regions, such as the hypothalamus, amygdala, and 
hippocampus, we conducted a preliminary follow-up func- 
tional magnetic resonance imaging (fMRI) study to explore 
whether rs6989250 is associated with neural activation in 
the stress-responsive brain circuits of cocaine-dependent 
subjects. 



Materials and methods 

Participants. Male and female subjects between the ages 
of 21 and 50 seeking inpatient treatment for cocaine 
dependence were evaluated for study participation (male: 
A/=36; female: A/=31). Subjects who met the Diagnostic 
and Statistical Manual of Mental Disorders (DSM-IV) criteria 
for current cocaine dependence were admitted to the CNRU 
(Clinical Neuroscience Research Unit) of the CMHC (Con- 
necticut Mental Health Center) for 2-4 weeks of inpatient 
treatment and research participation in a three-session 
laboratory experiment. Individuals who met the DSM-IV 
criteria for current dependence on additional psychoactive 
substances other than alcohol and nicotine were excluded. In 
addition, individuals on medications for medical or psychiatric 



diseases and those in need of alcohol detoxification were 
excluded from the study. The study procedures were 
approved by the Human Investigation Committee of the 
Yale University School of Medicine, and all subjects provided 
written informed consent. 

Procedures. During their first week of being an inpatient at 
the CNRU, the subjects were interviewed using the Struc- 
tured Clinical Interview for DSM-IV (SCID-I) 33 to determine 
psychiatric diagnoses. Baseline demographics and drug use 
history were also assessed. During the second week, the 
subjects participated in personalized script driven guided 
imagery procedures as outlined in previous work. 2 ' 3 ' 34-36 
In an imagery script development session, each subject 
identified a highly stressful event from their lives (rated by the 
subject as an 8 or above on a 10-point Likert scale for 
stressfulness), a personal drug cue-related event that involved 
people, a place and objects related to cocaine use, and a 
personal neutral relaxing event. Between weeks 3 and 4, the 
subjects participated in a habituation and imagery training 
session, followed by a 3-day laboratory experiment upon 
awakening from 0800 to 1030h, in which they were exposed 
to a stress, drug cue and neutral relaxing imagery provocation, 
one script per day, order randomized and counterbalanced 
across subjects, as described in our previous studies 3,37 
Script exposure occurred via audiotape and the length of each 
taped script was ~4.5 min, with a total imagery period of 5 min 
(see Supplementary Information for details). 

Laboratory measures. On each day of the experiment, 
cocaine craving was repeatedly assessed using a 10-point 
VAS (visual analog scale) that measured responses to the 
item 'desire for using cocaine at this moment'. Repeated 
blood samples were also collected for Cortisol assessments 
and processed at the YCCI (Yale Center for Clinical 
Investigation) Core Laboratories using standard radio- 
immunoassay procedures. On day 1 of the experiment, 
10 ml of blood was drawn for genotyping. 

Cocaine craving and plasma Cortisol responses were 
assessed at baseline, immediately following exposure to the 
imagery (time point 0) and at repeated recovery time points 
every 15 min (+15, +30, +45, +60 and +75) after 
exposure to imagery. All time points were utilized to assess 
response to brief stress, drug cue and neutral cue provocation 
in the laboratory. 

Prospective follow-up of cocaine relapse post-inpatient 
treatment. Appointments were scheduled for all participants 
for follow-up interviews in 14, 30 and 90 days after discharge 
from the inpatient unit. The patients were assessed 
for cocaine and other drug use with the Substance Use 
Calendar and urine and breathalyzer samples at each 
appointment, as described in our previous studies 3 ' 18 to 
prospectively assess cocaine relapse rates for the 90-day 
follow-up period (see Supplementary Information for follow- 
up assessments of relapse). 

OPRK1 genotyping. Genotyping was performed using the 
Addiction Array as described previously. 38 Briefly, haplotype 
tags for 1350 loci from 130 candidate genes were included in 
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this array. Additional 186 ancestry informative markers 
were included in the array to test population stratification. 
Population substructure was defined by applying the 186 
ancestry informative markers into the software 'Structure' 
(http://pritch.bsd.uchicago.edu/software.html). 39 Genotyping 
of rs6989250 was obtained from the Addiction Array based 
on lllumina GoldenGate Technology. The frequency of the 
G allele is 10% in African-Americans and 0% in Caucasians 
in our studied sample population. In all, 10 of the subjects 
in the study were CG heterozygous and 57 subjects were 
CC homozygous. 

Follow-up preliminary fMRI study. A follow-up preliminary 
fMRI study was also conducted in a demographically 
matched sample of five CG and eight CC patients (see 
Supplementary Methods for details). The magnetic reso- 
nance images were obtained using the 3-T Siemens Trio MRI 
system (Siemens Medical Solutions, Erlangen, Germany) 
with a single-channel, standard-quadrature head coil and a 
T2*-sensitive gradient-recalled single shot echo-planar pulse 
sequence. As in our previous fMRI studies, 36,40 a block 
design was used with six fMRI trials of personalized stress, 
drug cue and neutral scenarios (two separate events for each 
condition) that were presented in a randomized and counter- 
balanced order, and each script was presented via audio- 
tape. Each trial took 5min, which included a 1.5-min quiet 
baseline period followed by a 2.5-min imagery period (2 min 
of read imagery and 0.5 min of quiet imagery) and a 1-min 
quiet recovery period. During the baseline period, the 
participants were asked to remain in the scanner without 
any mental activity. Craving ratings were collected using 
a 10-point verbal Likert scale before and after each trial. 
To normalize any residual cravings from the previous trial, 
participants were engaged in 2 min of progressive relaxation 
between each trial (see Supplement Information). 

Statistics. The x 2 test and Student's Mest were applied to 
check for differences in the age, gender, years of education, 
duration and quantity of cocaine, alcohol and nicotine use 
between the two genotype groups (CG versus CC). 

Linear mixed effect models were employed for analysis of 
cocaine craving and Cortisol responses during the laboratory 
experiment using the SAS software package (Version 9, 
2006; SAS Institute, Cary, NC, USA). The between-subject 
factors for Group 2 (genotyping CG versus CC), the within- 
subject factors of Condition 3 (stress, drug cues and neutral) 
and Time point (baseline, 0, +15, +30, +45, +60 and 
+ 75) were the fixed effects, whereas Subject was the random 
effect factor in the model. Genetically inferred ancestry 
obtained from the 186 ancestry informative markers was a 
covariate in the model. Cohen's dwas used to estimate the 
effect size. The relationship between genotypes (CG versus 
CC) and relapse risk was examined using Log-rank survival 
analysis. 

fMRI analysis. Using XMedCon, 41 the raw imaging data 
were converted from the Digital Imaging and Communication 
in Medicine format into the Analyze format. The first 10 
images were removed from the beginning of each functional 
run to achieve a steady-state equilibrium between radio 



frequency pulsing and relaxation. The images were slice 
time-corrected and motion-corrected for three translational 
and three rotational directions; 42 any trial with a linear motion 
exceeding 1 .5 mm or a rotation of >2° was discarded. On an 
individual basis, a General Linear Model for each voxel in the 
brain volume was implemented with a regressor (time during 
imagery) for each trial per condition. In the General Linear 
Model, drift correction was included, such that drift regres- 
sors were used to remove the mean time course, linear trend, 
quadratic trend and cubic trend for each functional run. Due 
to the possibility of carryover effects from the imagery period, 
the recovery period was excluded from the data analysis. 
Each functional run was spatially smoothed using a 6-mm 
Gaussian kernel and normalized to generate beta maps 
(3.44 mm x 3.44 mm x 4 mm). 

Three registrations were sequentially implemented using 
the Yale Biol mage Suite software package (http://www. 
bioimagesuite.org) to adjust for individual anatomical differ- 
ences, 43 the linear registration of functional data into 2D 
anatomical images, the 2D to 3D (1 x 1 x 1 mm) linear regi- 
stration and non-linear registration to a reference 3D image. 
The reference image used was the Colin27 Brain 44 from the 
MNI (Montreal Neurological Institute) space. 

For the second-level group analysis, Analysis of Functional 
Neurol mages (AFNI) software 45 was utilized to run random 
mixed effect models. A 2 x 3 mixed effect model (group by 
condition) was conducted with Group as the between-subject 
fixed-effect factor; Condition (stress, drug cue and neutral) as 
the within-subject fixed-effect factor; and Subject as the 
random-effect factor. For further comparison, two-tailed 
Mests were conducted on the mean signal change contrasts 
(beta). To correct for multiple comparisons, a family-wise 
error (FWE) rate correction was applied via a Monte Carlo 
Simulation using Alpha Sim in AFNI. 46 To selectively isolate 
brain activity only specific to stress and drug components of 
imagery processing, the fMRI data results are presented for 
each contrast (stress-neutral and drug-neutral) using a within- 
subject subtraction method. 



Results 

Demographic data and individual characteristics by 
OPRK1 genotype. We analyzed 67 African-American 
cocaine-dependent individuals genotyped for rs6989250 
(CG = 10; CC = 57). There were no significant differences 
in age, years of education, the duration and amount of 
cocaine and alcohol use between the two genotype groups 
(CG versus CC) (P>0.05) (Table 1). However, the CG group 
included significantly more cigarette smokers than the CC 
group (100 versus 80%; P<0.05). Cigarette smoking was 
included as a covariate in the analyses. 

Cocaine craving. We first analyzed association of each 
SNP on the OPRK1 with cocaine craving. Among a total of 
11 SNPs of OPRK1 from the Addiction Array, the SNP, 
rs6989250, was the only SNP significantly associated with 
cocaine craving (P=0.02 uncorrected). There was no 
significant association of rs6989250 with craving after multi- 
ple test correction. However, considering that environmental 
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Table 1 Demographic characteristics of the two genotype groups 



Study 1: Laboratory experiment Study 2: fMRI 



Subject variable 


CG (N = 10) 


CC(N = 57) 


CG (N = 5) 


CC (N = 8) 


Gender (male) (%) 


5 (50%) 


31 (54.39%) 


2 (40%) 


3 (37.5%) 


Age (year±s.e.) 


37.7 (7.30) 


37.51 (6.48) 


38.2 (7.12) 


38.12(5.96) 


Average years of education ( ± s.e.) 


11.6 (0.97) 


12.4(1.91) 


11.8 (1.09) 


12.5 (1.60) 


Average years of cocaine use ( ± s.e.) 


1 1 .5 (5.97) 


10.05 (5.58) 


13(3.32) 


12.13(6.45) 


Average days of cocaine use/month ( ± s.e.) 


17.5(10.49) 


17.65 (8.10) 


15.8 (11.9) 


16(7.8) 


Total amount (g) of cocaine use in the past 90 days ( ± s.e.) 


29.23 (21.42) 


61.13 (106.2) 


20.44(15.65) 


56.17(42.08) 


Average years of alcohol use ( ± s.e.) 


15.89 (9.7) 


11.19 (8.34) 


14.75 (10.72) 


13.12(9.30) 


Average days of alcohol use/month ( ± s.e.) 


7.9 (7.4) 


10.55 (9.16) 


5.2 (6.3) 


9.25 (9.84) 


Regular cigarette smokers: N (%)* 


10(100%) 


46 (80%) 


5 (100%) 


6 (75%) 



Abbreviation: fMRI, functional magnetic resonance imaging. 

*P<0.05. Note: No significant differences were observed across all variables between the genotype groups in either the laboratory study or the fMRI study of the 
participants (Ps>0.05). The CG group included significantly more cigarette smokers than the CC group (P<0.05). 
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Figure 1 Mean (and standard error) drug craving response (averaged 
across all time points) during each of the neutral (N), stress (S) and drug cue (T) 
conditions for the two OPRK1 rs6989250 genotypes is shown. The genotype- 
by-condition two-way interaction was significant (P= 0.0056, Cohen's c/=0.79). 
Carriers of the CG genotype showed markedly stronger craving when exposed to 
the stress condition compared with individuals with the CC genotype (P= 0.0191, 
Cohen's c/=0.82). There was no significant difference detected between the two 
genotype groups in either the neutral (P= 0.0823) or drug cue conditions 
(P= 0.678). 
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Figure 2 Mean (and standard error) plasma Cortisol levels averaged across all 
time points in the neutral (N), stress (S) and drug cue (T) conditions for the OPRK1 
rs6989250 genotype groups are shown. There was a significant genotype-condition 
interaction (P=0.0003, Cohen's c/= 1 .01 ), with CG individuals showing higher 
Cortisol levels during stress compared with the neutral and drug cue conditions 
(P<0.0001 for stress versus neutral Cohen's cf=1.46; P= 0.003, Cohen's 
d= 1 .04 for cue versus neutral), whereas CC individuals showed no stress-induced 
increases in Cortisol compared with their Cortisol responses in the drug cue and 
neutral conditions (Ps>0.05). 



factor may manifest a marginal main gene effect, we focused 
on the interactive effect analysis of only rs6989250 with 
Condition. Consistent with our previous findings, there was a 
significant effect of Condition on drug craving (F= 20.04, 
df = 2, P< 0.001). Craving levels were higher under the 
stress and cue conditions compared with the neutral 
condition (P= 0.001). As mentioned above, there was no 
main effect of rs6989250 on drug craving after multiple test 
correction, but there was a significant genotype and condition 
interaction (F=5.40, df = 2, P= 0.006; effect size Cohen's 
c/=0.79). Individuals with a CG genotype exhibited a 
significantly stronger craving than subjects with a CC 
genotype during the stress condition (t= 2.370, df = 124, 
P=0.019, Cohen's d=0.81), but not in the neutral 
(P= 0.0823) and cue conditions (P= 0.678) (Figure 1). 



Plasma Cortisol. Because rs6989250 showed interactive 
effects with condition, we tested the interaction of rs6989250 
with plasma Cortisol level. As expected, we found a signi- 
ficant main effect of Condition on plasma Cortisol levels 
(F= 7.690, df = 2, P= 0.0007). The Cortisol levels were 
significantly higher during the stress and cue conditions 
compared with the neutral condition (f=3.79, df = 2, 
P= 0.0002 for stress versus neutral; f=2.74, df = 109, 
P= 0.0072 for cue versus neutral). Again, there was no 
significant main effect of genotype on Cortisol levels, but 
a significant genotype x condition interaction effect was 
observed (F=8.76, df = 2, P= 0.0003 Cohen's d=1.01). 
The CC group showed no significant increase in Cortisol 
during stress or drug cue exposure, while the CG group 
showed a highly significant increase in their Cortisol 
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response during stress (f=4.25, df = 109, P< 0.0001, Cohen's 
d= 1.46) and drug cue (f=3.04, df = 109, P= 0.003, Cohen's 
d=1.04) relative to the neutral condition (Figure 2). 

Cocaine relapse risk. We compared the subsequent 
cocaine relapse rates during the 90-day follow-up period 
after discharge from inpatient treatment and research for the 
CG and CC groups. The CG group showed a significantly 
greater relapse risk than the CC group (Log-rank x 2 = 7.157; 
df = 1 ; P= 0.0075) (Figure 3). For example, the abstinence 
rate in the first 14 days was only 10% in the CG group 
but 45% in the CC group. Thus, individuals with a CG 
genotype were more likely to relapse and return to cocaine 
use than the CC group, suggesting that the OPRK1 
rs6989250 SNP G allele may increase vulnerability to 
cocaine relapse. 

fMRI results. Whole-brain voxel-wise analysis indicated a 
significant genotype effect on the brain's response to stress 
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Figure 3 Log-rank survival analysis for abstinence versus cocaine relapse in 
the CG and CC genotype groups. The survival rate for the CG genotype group was 
significantly lower than for the CC genotype group (Log-rank % 2 = 5.242; 
P= 0.0075). Within 14days, only 10% individuals with the CG genotype remain 
abstinence of cocaine use; 45% individuals with the CC genotype remain 
abstinence. The /axis shows the proportion of participants who did not relapse; the 
Xaxis shows the 90-day follow-up time period. 



and to drug cues relative to neutral cue exposure in cocaine- 
dependent subjects (P<0.05, whole-brain FWE corrected; 
Table 2; Figure 4). In both the stress and the drug cues 
relative to neutral cue contrasts, the CG group showed 
greater activity in the hypothalamus, midbrain (ventral 
tegmental area, locus coeruleus), thalamus and cerebellum 
than the CC group. Additionally, the CG group displayed 
greater stress-induced activity in the right amygdala/hippo- 
campus and periaqueductal gray area (PAG) of the midbrain 
than the CC group (P<0.05 FWE corrected). 

Discussion 

This study showed that the putative functional SNP in the 
5'-regulatory region of OPRK1 is associated with significantly 
higher stress-induced craving, Cortisol responses and also 
increases subsequent risk of cocaine relapse in an African- 
American population with cocaine dependence. In a follow-up 
study with a subset of subjects, the CG genotype group 
exhibited stress- and drug cue-induced hyperactivity in 
hypothalamus/thalamus-midbrain-cerebellum regions, rela- 
tive to the CC group. Building on our previous findings 
identifying clinical and biological factors that are predictive of 
addiction relapse, 1 these results support the hypothesis of 
an interaction between an OPRK1 variant and stress reacti- 
vity that results in increased stress-induced cocaine craving, 
Cortisol responses and susceptibility to cocaine relapse 
in the African-American sample of cocaine-dependent 
individuals. 

Across different drug abusing samples, including cocaine, 
previous research has shown that exposure to acute stress 
promotes drug craving and significantly increases the risk of 
early relapse. 47-49 Current results indicate that individual 
variation differed by the OPfl/Ogenotype significantly con- 
tributed to the individual variation in the stress-induced 
craving responses. During laboratory-induced stress expo- 
sure, subjective craving was remarkably higher in CG 
subjects than in CC subjects. In previous research, we have 
shown that stress-induced drug craving and stress-induced 
Cortisol responses were predictive of earlier relapse to 
cocaine use, 3 and thus we expected that such higher stress 
and drug craving responses in the CG group would make 
them more susceptible to relapse. Indeed, we find in the 
current sample that 90% of the CG subjects relapsed within 
the first 2 weeks after treatment, suggesting that these 



Table 2 Brain regions showing group differences between the CG group and the CC group (CG>CC) 



Stress — neutral 



Coordinates Volume 



Regions of activation 


Lat 


X 


Y 


Z 


(mm 3 ) 


Amygdala/hippocampus 


R 


25 


-7 


-18 


1226 


Hypothalamus/thalamus 


B 


-1 


-5 


-7 


1735 


Midbrain/brain stem 


B 


1 


-26 


-27 


6925 


Cerebellum 


B 


0 


-39 


-28 


4823 



Drug cue — neutral 
Coordinates Volume 



t 


d 


X 


Y 


Z 


(mm 3 ) 


t 


d 


2.31 


1.32 














2.39 


1.36 


5 


-8 


-6 


738 


2.2 


1.25 


2.38 


1.36 


1 


-28 


-27 


6832 


2.5 


1.43 


2.46 


1.40 


2 


-43 


-27 


6409 


2.5 


1.43 



Abbreviations: B, bilateral; d, Cohen's d post hoc estimates of effect size; Lat, laterality; R, right. 
Note: Significant activations at P<0.05 (two-tailed, whole-brain FWE corrected). 
Montreal Neurological Institute (MNI) coordinates were used. 
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Figure 4 Whole brain analysis of neural activation during (a) stress and (b) drug cue exposure in the CG group relative to the CC group. During the stress and drug cue 
relative to the neutral condition, the CG group showed increased activity in the hypothalamus, midbrain (ventral tegmental area and locus coeruleus), thalamus and cerebellum 
relative to the CC group. Additionally, the CG group displayed increased stress-induced activity in the right amygdala/hippocampus and periaqueductal gray (PAG) area of the 
midbrain compared with the CC group (P<0.05, whole-brain family-wise error (FWE) corrected). The blood oxygenation level-dependent (BOLD) functional magnetic 
resonance imaging (fMRI) signal in each significant cluster (Table 2) is displayed in bar graphs for each genotype group in (a) stress-neutral and (b) drug-neutral conditions 
(see Table 2 for f-values and Cohen's d). Coordinates are given in Montreal Neurological Institute (MNI) space. 



individuals, who showed greater stress-induced cocaine and was not seen in response to drug cues. These data 
craving, are more vulnerable to an earlier risk of relapse. suggest that the OPRK1 variant specifically affects stress 
Interestingly, this effect was specific to the stress condition sensitization of drug craving and increases the risk of relapse 
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for individuals with this genetic susceptibility. Although these 
findings were generated from a small sample size, we see 
medium effect size in the Genotype x Condition interaction 
effect for cocaine craving responses. 

Consistent with the evidence that KOR modulates stress via 
the hypothalamic-pituitary-adrenal axis, our results demon- 
strate that a genetic variant of OPRK1 influences Cortisol 
levels during stress. 50 ' 51 Previous studies show that KOR 
agonists increase corticosterone in rats, 52 elevate plasma 
Cortisol levels in humans 53 and stimulate plasma ACTH and 
Cortisol release in a dose-dependent manner in non-human 
primates 54 The stimulatory effects of the KOR agonist on the 
hypothalamic-pituitary-adrenal axis were blocked by the 
KOR antagonist nor-binaltorphimine. 54 However, the endo- 
genous kappa opioid peptide dynorphin increases the release 
of corticotrophin-releasing factor in the hypothalamus. 55 More 
importantly, the results of the present study indicate that 
genetic variation in OPRK1 is associated with individual 
hypothalamic-pituitary-adrenal responses to stress, as 
shown by the effect of rs6989250 on the Cortisol response in 
the studied population. Here, we found an interactive effect on 
Cortisol of the rs6989250 genotype and the experimental 
provocation of stress and drug cue conditions. The Cortisol 
level was significantly higher under stress and drug cue 
conditions relative to the neutral condition among the subjects 
with the CG genotype, an effect similar to that seen with a 
KOR agonist. The genotype effect on Cortisol response was 
large in both stress and drug cue conditions, which is reflected 
by Cohen's d>1. Although the function of rs6989250 is 
unknown, the G allele creates an NF-kB transcription factor 
binding site, presumably resulting in higher OPRK1 expres- 
sion and increased hypothalamic-pituitary-adrenal axis 
function. Thus, the current findings suggest that individuals 
with the G allele are more vulnerable to stress and stress- 
induced Cortisol response. 

The two genotype groups showed significant differences in 
limbic-midbrain regions, specifically in hypothalamus/ 
thalamus-midbrain-cerebellum regions during stress and 
drug cue exposure. The CG group showed greater activity in 
the hypothalamus, midbrain, thalamus and cerebellum than 
the CC group during the stress and drug cue relative to the 
neutral cue trials. KORs are primarily expressed in the 
hypothalamus and midbrain regions, 56 consistent with our 
findings of hyperresponsivity to negative emotion/stress and 
drug cue challenge in the CG group possibly reflecting an 
upregulated neuronal function in the regions that are densely 
packed with KORs. Specifically, the midbrain VTA is important 
in reward and incentive processing via its involvement in 
dopamine firing and KOR control of the VTA dopamine 
projection to the prefrontal cortex. 57 Additionally, the CG 
group showed stress-induced hyperactivity in the right 
amygdala/hippocampus and midbrain PAG. The amygdala 
and hippocampus have been associated with emotional and 
stress processing, and hyperactivity in these regions has 
been reported during exposure to acute and chronic 
stress. 46 ' 58 Specifically, right-lateralized limbic activity has 
been observed during aversive emotional processing, such as 
in response to fear and negative emotion. 59,60 The PAG, a 
midbrain region with abundant connections to the amygdala, 
is involved in pain, fear and anxiety, 61 and PAG opioid 



receptors modulate aversive emotional processing. 6 
The stress-induced hyperactivation observed in the right 
hypothalamus/thalamus regions and the PAG in the CG 
group is consistent with previous findings of association 
with kappa opioid gene expression related to stress sensitivity 
and dysphoria. 12 These results suggest that the neuronal 
function in limbic-midbrain regions, which are key areas 
that contain densely distributed KORs, are upregulated 
during stress and reward processing in CG relative to 
CC individuals, which could make them more vulnerable to 
stress and drug craving responses as seen in the laboratory 
experiment. 

It should be noted that genotype effects were of moderate to 
large effect sizes for each of the dependent measures in this 
study and larger than those reported for common variants in 
complex traits. This large genotype effect may be attributable 
to careful endophenotype characterization of the stress 
response and minimal population admixture that potentially 
reduced sample heterogeneity. It is also possible that 
rs6989250 is indeed functional, or is in linkage disequilibrium 
with a functional locus on OPRK1. 

A limitation of the current findings is the small size of the 
study sample. Also, the functionality of rs6989250 has not 
been evaluated. We were unable to perform haplotype 
analysis because of insufficient statistical power. With these 
limitations in mind, we suggest caution in interpreting current 
results. A larger sample size from African-Americans and from 
other ethnicities and genotyping higher SNP density will be 
needed in future studies to replicate the results of this study. 
Further evaluation of the function of the 5' regulatory variant 
will also provide more insight into the mechanisms underlying 
the role of OPRK1 in modulating stress, drug craving and 
relapse risk associated with cocaine dependence. For 
example, deep sequencing of the 5' OPRK1 variant in 
cocaine-dependent individuals may identify uncommon or 
rare causal variant(s) of OPRK1. 

In summary, to our knowledge, this is the first report 
of a putative functional SNP in OPRK1 that is linked to 
individual variation in stress-related drug craving, Cortisol and 
cocaine relapse risk as well as hyperactivation of hypothala- 
mus/thalamus regions during stress and drug cue exposure in 
cocaine-dependent subjects. Our results indicate a role for 
OPRK1 in increasing stress, drug craving and cocaine relapse 
risk in African-Americans. These results also support the need 
for individualized treatment of cocaine dependence, wherein 
assessment of genotype and stress vulnerability may help 
identify those individuals who are most susceptible to early 
relapse so that development and testing of new, specialized 
treatments may be considered for such individuals. 
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